Similar glassy features in the NMR response of pure and disordered Lai gsSro.^CuC^ 
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High T c superconductivity in La2_ I Sr a! Cu04 coexists with (striped and glassy) magnetic order. 
Here, we report NMR measurements of the 139 La spin-lattice relaxation, which displays a stretched- 
exponential time dependence, in both pure and disordered x = 0.12 single crystals. An analysis in 
terms of a distribution of relaxation rates 139 j 1 l _1 indicates that i) the spin-freezing temperature is 
spatially inhomogeneous with an onset at y° nsct = 20 K for the pristine samples, and ii) the width 
of the T] -1 distribution in the vicinity of y° nsct ; s insensitive to an ~1 % level of atomic disorder 
in CuC>2 planes. This suggests that the stretched-exponential 139 La relaxation, considered as a 
manifestation of the systems glassiness, may not arise from quenched disorder. 



I. INTRODUCTION 

The coexistence of magnetic order with superconduc- 
tivity is a prominent feature of L^-^Sr^CuO^ [T] and 
other underdoped high T c cuprates [3], even in zero 
magnetic field. However, neither the origin of this static 
magnetism nor the reason for its glassy character are 
fully understood. An important school of thought fo- 
cuses on stripe physics jlj. Indeed, the temperature of 
magnetic freezing, T g , in La2- 2; Sr 2 ;Cu04 is peaked in 
the vicinity of x ~ | (Fig. [TJ, and neutron scattering 
studies [5J |7] reveal long-range antiferromagnetic or- 
der with the same typical modulation as in the materials 
presenting direct evidence for charge stripe order [8] . An- 
other approach relies on electronic and magnetic inhomo- 
geneities generated by quenched disorder [5] . Undoped or 
weakly hole-doped droplets may form the magnetic clus- 
ters which freeze at low temperature, as suggested by 
the recent nuclear magnetic resonance (NMR) evidence 
for a nanoscale inhomogcneity of the hole concentration 
in La2- a; Sr a; Cu04 Furthermore, the importance of 

quenched disorder could be supported by the glassy fea- 
tures observed in superconducting samples (x > 0.06) [1J. 
These are reminiscent of the spin-glass behaviour, well 
documented for 0.03 < x < 0.05. Spatial heterogeneity 
is particularly evident in a number of magnetic measure- 
ments, such as the stretched-exponential NMR relaxation 
of 139 La nuclei. This heterogeneity develops as the mag- 
netic fluctuations slow down over a substantial temper- 
ature range on cooling above the freezing temperature 
T g , the value of which depends on the timescale of the 
measuring probe. These properties are typical of glassy 
systems. 

Glassiness is, however, also present in materials where 
stripe order is well-defined and relatively long range [5]. 
Furthermore, from the T dependence of an average 139 La 
spin-lattice relaxation rate T^ 1 , Curro et ai. [12] con- 
cluded that the distribution of magnetic correlation times 
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FIG. 1: (Color online) Sketch of the phase diagram of 
La2-iSr a; Cu04. AF denotes antiferromagnetic, SC supercon- 
ducting, and SG spin glass phase. T g defines the temperature 
of the transition to the frozen state which has been named a 
"cluster spin glass" before the discovery of stripes (see [T] and 
Refs. therein). T g is defined from either NMR or /iSR mea- 
surements. These two low energy magnetic probes of compa- 
rable (though not necessarily identical) timescales give sim- 
ilar values of T g , enabling one to obtain a coherent phase 
diagram [TJ. Remarkably, T g (together with the spin stiff- 
ness - see text) is enhanced for x ~ 0.12. One of the ques- 
tions regarding the cluster spin-glass phase is whether the 
non-monotonic behavior of T g vs. x suggests different physics 
for x ~ 0.12 than for x < 0.10. Note that the main difference 
between L^-zSr^CuO^ and compounds presenting direct ev- 
idence for charge stripe order (e.g. La2- ;c Ba I Cu04) is that 
the peak of T g vs. x around x ~ 0.12 is considerably wider 
in these latter compounds [3]. 



is similar in several materials with very different hole or 
impurity contents. This has led to the suggestion that 
glassiness in these materials is not due to quenched disor- 
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der but is self-generated by the charge stripes [13] ■ Later, 
Hunt et al. |14j determined the distribution of correlation 
times in stripe-ordered materials in a more direct way 
than Curro et a?.[T2], i.e. from the recovery of the 139 La 
NMR signal. The stretched-exponential relaxation was 
thus better characterized, but its origin was not the fo- 
cus of the discussion. 

In this paper, we report 139 La NMR measurements in 
Lai.88Sro.i2Cu04, where T g is maximum, as illustrated 
in Fig. [T] We quantify the 139 La stretched-exponential 
relaxation and discuss the issues of inhomogeneity and 
atomic disorder. 
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II. SAMPLES 

Single crystals, of typical size 2x3x4 mm 3 , were cut 
from rods grown by the traveling solvent floating zone 
method. Growth at a rate of ~1 mm/h yielded "pristine" 
samples with the standard T c of 30 K [5], while growth 
at 0.2 mm/h produced a "disordered" sample with a low 
T c of 10 K (T c onsot = 12 K) [15]. Neutron scattering (NS) 
studies have shown that the two crystals have the same 
structure as well as an identical temperature of the tran- 
sition from the high T tetragonal (HTT) phase to the low 
T orthorhombic (LTO) phase. These studies have also re- 
vealed an identical incommensurability of the magnetic 
peaks. Since these quantities depend strongly on x, the 
results demonstrate that the two crystals have the same 
doping x = 0.12. On the other hand, magnetic Bragg 
peaks appear below Tf s = 30 K [5] and 25 K [H] for 
the pristine and disordered samples, respectively. All of 
these properties, as well as an upturn of the zero-field 
in- plane resistivity below 80 K, are consistent with the 
presence of ~ 1 % of non-magnetic defects in the disor- 
dered sample [T5] . These defects likely correspond to Cu 
vacancies, which produce the same magnetic effects as 
non-magnetic impurities |17[ 118] . 



III. NMR METHODS 

As in [T5], the applied magnetic field was tilted away 
from the c-axis by 9 ~ 10°, in order to obtain a sharp 
peak on the high frequency edge of the 139 La central, 
(+\ <-> — \), Zccman transition. T\ was measured on 
this peak, shown as a shaded part of the sample spectra 
in the inset to Fig. [2] The T dependence of T\ is identical 
on other peaks in the central transition spectrum. Ex- 
periments were performed in fields of 9 T (for x = 0.12, 
T c = 10 K sample) and 14 T (for x = 0.12, T c = 30 K, 
and for x — 0.10 samples). 

In Fig. [2] we plot the temperature dependence of 139 La 
linewidth for the high frequency peak (shaded region in 
the inset) of (+5 — h) transition for Lai.8sSro.i2Cu04 
(T c = 30 K sample). The linewidth broadens with de- 
creasing temperature indicating increasing distribution 



FIG. 2: (Color online) The temperature dependence of 139 La 
linewidth (filled symbols) for the high frequency peak (shaded 
region in the inset) of (+f <_> — |) transition in the mag- 
netic field applied at ~ 10° away from the c-axis for 
Lai.88Sr .i2CuO4 (T c = 30 K sample). The solid line is 
guide to the eye. Inset: Normalized spectra at T = 100 K 
and T = 7.1 K as denoted. The shaded areas mark the 
high frequency peak where the rates presented here are de- 
termined. The low frequency peak contains the signal from 
other quadrupolarly split (+~ «-» — |) transition lines |19| . 

of local static magnetization as the temperature is low- 
ered. 

We remark that even at the lowest T the signal from 
other central transition lines is insignificant at the fre- 
quency of the the high frequency peak and its high fre- 
quency side, for this orientation of the applied field. 
Therefore, the temperature evolution of the measured re- 
laxation rates is intrinsic to magnetism and not a result 
of signal overlap from other central transition lines. 

We also point out that no sign of phase separation 
is detected in our NMR data, in apparent disagreement 
with /iSR data showing a magnetic volume fraction of 
~20% for x = 0.12 [20]. Since NMR was performed here 
in a high magnetic field, part of the discrepancy might be 
resolved by an increase of the magnetic volume fraction 
with the applied field [2T]. |2"2] . 

IV. ANALYSIS OF Ti DATA 

In Figure [3] we display the time (t) dependence of the 
139 La longitudinal magnetization A4(t) = M z (t)/M z (oo) 
after a comb of ^ saturation pulses, in the pristine sam- 
ple. Clearly, the data cannot be fitted by the theoretical 
formula |23j , which corresponds to a = 1 ( "exponential" 
relaxation) in 

Mo^Tf 1 ) = 1 - 0.714 e~( 28 ^) - 0.206 e~( 15 ^) 

-0.068e~( 6 ^) - 0.012e _ (^) . (1) 

We point out that the above statement is true for data 
at any temperature below w 80 K, as evident in Fig. |3j; 
when a < 1. In La2-xSr a: Cu04 materials, the stretched 
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FIG. 3: (Color online) (a) An example of the time de- 
pendence of the 139 La normalized longitudinal magnetization 
M(t) = M z (t)/Mz(oc) after saturation of the central tran- 
sition (dots) in Lai.88Sro.i2CuC>4 (T c = 30 K sample). The 
thin black line is a fit to Eq. 1 with a = 1. The blue (dark 
gray) dashed line is a fit to Eq. 1 with a — 0.48. The thick 
red (light gray) line is the fit to Eq. [2] (b) and (c) Compar- 
isons of results obtained from the stretched-exponential fit 
(open symbols) and from the Gaussian distribution fit (filled 
symbols). See text for details. 



relaxation, i.e. the deviation from A4 a =i{t, T^ 1 ), has 
been attributed to a distribution of 139 T 1 _1 values. We 
shall comment on this interpretation below. 

In an ideal case, one-to-one correspondence between 
dynamic and static inhomogcneities can be revealed by 
measuring Tj" as a function of position across the NMR 
line shape (see [TTJ [51] for example). However, despite 
the observation of a continuous line-broadening (Shown 
in Fig. 2), no significant frequency dependence of T-J -1 
was found across the 139 La line. Thus, to quantify the 
inhomogeneities one must resort to alternative analysis 
of M.(t) data, in two possible ways: 

(i) First, a fit to Eq. [T]is made with the stretching ex- 
ponent ct^ 1. This provides a phenomenological account 
of the distribution of T\ values (Fig.|3|. 

(ii) Second, the formula Ai a =i(t, T± ), is convoluted 
with a chosen probability distribution function of T± . 
We found that a good fit to the data is provided by the 
Gaussian distribution on a logarithmic (logT^ 1 ) scale: 

M Q (t) = (7^72 Tiog)- 1 x 

-2(iogi? 1 -io g Tr 1 ) 2 /^o g M a=1 {t,Rx) d{\ogR{). (2) 



This fit is defined by only two parameters: the most prob- 
able relaxation rate Tj~ , (i.e., the center of the Gaus- 
sian), and the width of the distribution ai og on a logio 
scale. Using the logTj" 1 space is more physical when 
very broad distributions of relaxation rates are expected. 
This also naturally avoids the introduction of an artifi- 
cial low T^ 1 cut-off needed to eliminate unphysical neg- 
ative values encountered when the linear scale is used. 
The fit yields T-f 1 values equivalent to the stretched fit 
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FIG. 4: (Color online) The Gaussian probability distribution 
of the relaxation rate as a function of T for Lai.88Sro.i2Cu04 
(T c — 30 K sample). The rates are in units of s~ . The 
amplitude of distributions are normalized to one and curves 
offset for clarity. The vertical line indicates a typical maxi- 
mum value (T 1 ' 1 ) max ~ 10 3 s _1 in the system. Green color 
(lightest gray) depicts the T interval (from 100 K to 25 K) 
of slowing down, red (lighter gray) represents the T range 
(from 20 K to 9 K) where spins freeze throughout the sample, 
and navy blue (dark gray) is used when all of the sample is 
magnetically frozen on the NMR timescale. 



M a (t, Tj -1 ), as confirmed from the analysis of our data 
shown in Fig. [3|). Furthermore, the value of <7i og directly 
shows over how many orders of magnitude spreads the 
distribution of relaxation rates. Direct insight into this 
parameter is the main advantage of this fit (see Fig. EJ. 
The a\ og also appears to be linearly related to the value 
of a, as demonstrated in Fig. R5p, in agreement with pre- 
dictions of the recent theoretical work of Johnston for 
a > 0.5 [35]. We also remark that we cannot experi- 
mentally determine the unique/exact shape of the dis- 
tribution function. Equally good fits of the data can be 
achieved by assuming, e.g. an asymmetric Gaussian or 
a Lorentzian distribution for T-j - . Thus, we choose to 
analyze the data assuming a distribution function of the 
simplest form, a Gaussian. However, the results to be 
discussed below have been found to be insensitive to the 
choice of the exact form of the distribution. 



V. TEMPERATURE DEPENDENCE OF THE 
SPIN DYNAMICS 

In this section, we shall first give a qualitative descrip- 
tion of the temperature dependence of the mean relax- 
ation rate , in order to define the freezing tempera- 
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FIG. 5: (Color online) (a) 139 La spin-lattice relaxation 
rate Tj -1 defined as the center of the Gaussian distribution 
V, for the two x = 0.12 samples and for an x = 0.10 
sample, (b) P max = P(log(r- 1 ) max )/P(log(T 1 - 1 )) where 
(T 1 r 1 )max = 10 3 s _1 (same symbol code as in main panel). 
Vertical dashed lines denote T = 15 and 20 K, where the 
first spins become frozen in the samples, (c) Width <j\ oe of 
VQogTi) (same symbol code as in upper panel). The lines 
are guides to the eye. The downturn of o"i og at low T for 
x = 0.10 might be due to a loss of NMR signal. 



ture T g on the timescale of NMR. Then, tentative fits 
of the T dependence of T^ 1 will be described. This 
parametrization of the data enables primarily a com- 
parison between different La2-xSr x Cu04 samples, while 
it cannot providing a physical picture of the spin-glass 
freezing in the system. Furthermore, the shortcomings 
of such analysis based on the mean value of the distribu- 
tion of T-f 1 values are alluded to in the next section. 

The temperature dependence of T± may be understood, 
at least qualitatively, from the following expression: 



1 



2t c 



(3) 
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where r c is the correlation time, h± = (h%, 
component of the hyperfinc field perpendicular to the 
field direction, and u n the NMR frequency [55]. At high 
temperatures the correlation time is short, that is the 
condition Tf 1 ^> uj n is satisfied. As the dynamics of the 
system slows down on cooling, r" 1 decreases causing an 
increase of Tf 1 . This occurs down to the temperature 
T = T^ MR at which the condition rf 1 = cu n is reached. 

Upon further cooling, t~ x continues to decrease but T-f 
decreases. Thus, T^ M , the temperature of freezing on 
the timescale of 139 La NMR (w n ~ 10 8 Hz), is defined 
as the temperature at which the relaxation rate is at its 
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For our pristine x = 0.12 sample, the peak of the mean 

Tf 1 occurs at T g NMR = 13 K. The temperature at which 

the increase of T£~ becomes noticeable may be defined 

as T slow = 45 K. Interestingly, the ratio T slow /T g NMR ~ 

3.5 ± 0.5 is much lower for x = 0.12 than for x = 0.10 
( T s1ow/ t nmr > 7A j and for other valueg of x [ioJ [3TJ. 

This statement holds even if an onset T^ MR is considered 
as described in the next section. Note also that the ratio 
is believed to be magnetic field independent for fields 
(H sS 14 T) investigated here [H]. 

A more quantitative approach requires an analysis of 
the temperature dependence of T\. Although it is not 
clear on which theoretical model such an analysis should 
be based, two models which have been used in the context 
of spin-freezing in the cuprates can be used. 

First, we use Eq. [3] with an activated correlation time 
t c = to exp(E a /kBT) to fit the temperature dependence 
of Xf 1 data with To and E a as fitting parameters, as 
depicted in Fig. [6] This fit allows to extract the effec- 
tive 'energy barrier' (E a ) for the activation of a ther- 
mally driven spin freezing process. Fitting the data for 
t nmr < T < ysiow^ we find the effective energy bar- 
rier E a — 140 ± 30 K for our pristine x — 0.12 sample. 
This value is to be compared with E a — 84 ± 20 K for 
the disordered x — 12% sample and E a — 13 ± 3 K for 
x = 10%. 

Another possibility is to use the renormalized classical 
form of the 2D Heisenberg model on a square lattice, 
which also captures well the increase of T^ 1 with decreas- 
ing temperature above T'^ ,MR (see also [H]). This form 
allows to extract the effective spin stiffness (p s ) param- 
eter. As illustrated in Fig. [6] in the temperature range 
t nmr <t < T slow , the data for both x = 0.12 samples 
fit well to the low temperature limit (T < 2irp s /2) of the 
renormalized classical form, expressed as 



1 e c ' T ■ (T/Cfl 2 
T\ * C-(l + T/C) 3 ' 



(•5) 



Here C = 2-xp s is a fitting parameter. We find a value 
of p s — 25 ± 5 K for the spin stiffness of our pris- 
tine x = 0.12 sample. For the disordered x — 0.12 
sample, p s — 19 ± 4 K is lower. An even weaker value 
(p s = 2 ± 1 K) is found for the x = 0.10 sample. This 
correlates with the fact that T slow /T g NMR > 7.4 for this 
sample. Due to its weak spin stiffness parameter, the 
data for the x — 0.10 sample are fitted to the high tem- 
perature limit (T > 2Trp s /2) of the renormalized classical 
form [28], given by 
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with C = 27rp s as a fitting parameter. Not surprisingly 
(given the similarities of the fitting formulas in this tern- 



5 




20 40 60 80 100 

Temperature [K ] 



FIG. 6: (Color online) Fits to the relaxation rate T~ x data as 
plotted in Fig. [5^. The dashed lines are fits to the Eq. |3]with 
t~c ~ to exp(_E a /fesT). The solid lines are fits to the Eq. [5] 
for both x — 0.12 samples and to the Eq. [6j for x — 0.10, as 
described in the text. 

perature range), the extracted values of the spin stiff- 
ness 2irp s are comparable to the values of the effective 
energy barrier in all the three samples. On the other 
hand, the spin-stiffness values are smaller than those in 
La-based cuprates with LTT structure and charge-stripe 

order QUH5]. 

In conclusion, it appears that x = 0.12 is the concen- 
tration for which the temperature range of slowing down 
of the magnetic fluctuations is the least extended and/or 
the spin-stiffness is the strongest. We remark that the 
last point is in agreement with results in Ref. [25] . 



VI. INHOMOGENEITY OF T g 

For a more adequate description of the spin freezing, 
it is in fact the maximum, and not the mean, Tj -1 value 
present in the system which should be considered. We 
define this value to be the rate at which the probabil- 
ity distribution is 1% of its maximum. Based on re- 
sults shown in Figs. [4] and [5]), we deduce that a typical 
maximum value is approximately (Tj - ) max ~ 10 3 s _1 . 
It is apparent in Fig. [5b that the weight of the Gaus- 
sian distribution at T x = 10 3 s _1 becomes sizeable at 
T™ R (onset) = 20 K. This means that spins begin to 
be frozen well above 13 K in some regions of the sample. 
Our T^^ ct = 20 K agrees perfectly with the appearance 
of a coherent precession of the muon spin in a x = 0.12 
single crystal similar to ours [2D] . 

Here, it must be recognized that (T-f ) raax = 10 3 s _1 
probably also represents a cutoff value above which 
the 139 La NMR signal is unobservable (wipeout phe- 
nomenon) [T2] H3J [TpJ . However, the consistency of the 
above analysis [3D] suggests that our data correctly de- 
scribe the spin dynamics of the system at least down to 
20 K. 

For the disordered sample, the freezing occurs at a 



lower temperature: T 1 1 vs. T shows a peak at T — 10 K 
and T g NMR ' onsot = 15 K. This is consistent with the quali- 
tative analysis discussed in the previous section and with 
the decrease of T™ s from 30 K (pristine) [S] [BJ [7J to 25 K 
(disordered) [15] , as determined by neutron scattering. 
Since the elastic NS signal is integrated over a typical 
energy window of ~1 meV and the slowing down of the 
fluctuations occurs on a relatively wide T range, we ob- 
serve > T™f et . 



VII. INHOMOGENEITY OF THE SPIN 
DYNAMICS 

The distribution of T± values acquires a significant 
width only at T ~ 60 K and below. The fact that 
CTi g ~ at T = 80 K and 100 K is remarkable. This 
means that the distribution of T^ 1 values seen by 63 Cu 
or 17 NMR at these temperatures, and attributed to a 
nanoscalc inhomogeneity of the hole concentration [llj . 
is completely absent in the 139 La T\ data. This con- 
trasting situation is explained by two facts: First, the 
hyperfine field at 139 La sites results from the coupling 
to several electronic sites in two different Cu02 planes. 
Thus, the hyperfine field is spatially more homogeneous 
than the electronic density. Second, the relatively weak 
amplitude of this (transferred) hyperfine field produces 
much weaker differences in the T\ values for 139 La than 
for 63 Cu or 17 0. Still, the 139 Ti inhomogeneity, which 
shows up at lower temperature, appears to be correlated 
to inhomogcncitics in the static local magnetization, as 
revealed by the similar temperature dependence of both 
the distribution width <7i og and the static 139 La linewidth 
shown in Fig. [3J; and Fig. [2] respectively. However, for 
the reasons given above, no direct correlation between 
the inhomogeneity of 139 T -1 and the resonance frequency 
across the 139 La line could be established. 

Since 17 and 63 ' 65 Cu NMR signals undergo a signifi- 
cant wipeout below ^50 K [19 , 32J, it is difficult to know 
whether the inhomogeneity probed in 139 La T\ measure- 
ments below ^60 K is distinct from the one probed at 
higher temperatures and attributed to spatial variation of 
the hole concentration. At any given temperature below 
^20 K, the difference in ^T^ 1 values between x = 0.12 
and x — 0.10 samples exceeds an order of magnitude. 
However, even if the amplitude of the nanoscale doping 
inhomogeneity does not produce a significant effects on 
the 139 La NMR at ~100 K, the possibility that it pro- 
duces the large 139 Ti distribution at low T cannot be dis- 
carded. On the other hand, it is possible that magnetic 
heterogeneity (ubiquitous in glassy systems) develops 
prior to the glass transition in L^-zSr^CuO^ in addi- 
tion to a nanoscale electronic inhomogeneity. In support 
to this is the fact that spin freezing and substantial mag- 
netic inhomogeneity are reported in YBa 2 Cu 3 6+:i; [2J, 
while a significant nanoscale variation of the hole con- 
centration is apparently absent in this system 33J . 

As evident in Fig. [JjJd, T\ values are typically dis- 
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tributed over more than one order of magnitude at low T, 
that is eriog reaches a value of - 1.5 at T = T ff NMR . This 
width does not depend strongly on the number of holes 
since similar values of <r\ og are observed for x — 0.10 [TS] 
(Fig. |J and x = 0.06 [3T] (not shown) at T ~ T g NMR . 

Strikingly, the same value of cri og ~ 1.5 is found at 
T = T^mr in both pristine and disordered samples for 
x = 0.12 (Fig. [5t). Thus, the distribution width appears 
to be insensitive to an ~ 1% of in-plane disorder as well. 
One could argue that this result can be explained by the 
fact that Lai.ggSr .i2CuO4 is already a significantly dis- 
ordered material in its pristine version. However, this ar- 
gument docs not hold since a typical 1 % of non-magnetic 
impurities or vacancies is highly detrimental to T c and it 
clearly affects other magnetic properties as well [T6"| . 

Finally, we will comment on alternative explanations 
of the stretched-exponential behavior of the NMR relax- 
ation. It was suggested that the inhomogeneous mag- 
netic state is related to extended charge density waves 
with imaginary order parameter (id-CDW) [33]. Within 
this model, the observed inhomogeneity in the NMR and 
/iSR quantities would originate from sliding motions of 
orbital currents coexisting with the <i-wave superconduct- 
ing state. It is known that any order parameter of d-wave 
symmetry should be highly sensitive to impurities. This 
is in contrast to our data, showing that the inhomogene- 
ity of 139 Ti is insensitive to 1% of atomic disorder. There- 
fore, the id-CDW scenario is unlikely to account for our 
observations. 

Furthermore, a power-law time dependence of the spin- 
spin correlation function, instead of a distribution, is 
also debated in canonical spin-glasses as an explana- 
tion for the stretched-exponential NMR or fiSR relax- 
ation |35j . Anyhow, we did not consider this possibility as 
the presence of magnetic heterogeneity in L^-^Sr^CuC^ 



is established by various experimental facts, such as the 
partial wipeout of the NMR signal or the T-dependent 
broadening of the NMR lincwidth. 



VIII. CONCLUSION 

We have presented an NMR investigation of 
Lai.88Sr .i2CuO4, a prototypal material for study- 
ing spin-glass and stripe-ordering instabilities in the 
superconducting regime. Our analysis can be viewed 
as a parametrization of the stretched-exponential spin- 
lattice relaxation of 139 La nuclei. We observe that 
this phenomenon is not affected by a ~1 % level of 
disorder. This result might thus support proposals [T3] 
of heterogeneous dynamics, or more generally of glassi- 
ness, which is not due to quenched disorder. However, 
how much of the magnetic heterogeneity/glassiness of 
La2- 2: Sr 2 ,Cu04 is attributable to nanoscale variation 
of the hole doping remains unclear. It would also be 
interesting to investigate whether our results for the 
special x = 0.12 material still hold in samples with lower 
Sr concentration, including the non superconducting 
cluster spin- glass region of the phase diagram. The 
glassy nature of magnetism in superconducting cuprates 
clearly calls for further theoretical and experimental 
consideration, especially given the recent observation of 
an electronic glass by scanning tunneling microscopy [36J . 
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